\(\int \csc (e+f x) \sqrt {a+a \sin (e+f x)} \sqrt {c+d \sin (e+f x)} \, dx\) [35]

Optimal result

Integrand size = 35, antiderivative size = 123 \[ \int \csc (e+f x) \sqrt {a+a \sin (e+f x)} \sqrt {c+d \sin (e+f x)} \, dx=-\frac {2 \sqrt {a} \sqrt {d} \arctan \left (\frac {\sqrt {a} \sqrt {d} \cos (e+f x)}{\sqrt {a+a \sin (e+f x)} \sqrt {c+d \sin (e+f x)}}\right )}{f}-\frac {2 \sqrt {a} \sqrt {c} \text {arctanh}\left (\frac {\sqrt {a} \sqrt {c} \cos (e+f x)}{\sqrt {a+a \sin (e+f x)} \sqrt {c+d \sin (e+f x)}}\right )}{f} \]

[Out]

-2*arctanh(cos(f*x+e)*a^(1/2)*c^(1/2)/(a+a*sin(f*x+e))^(1/2)/(c+d*sin(f*x+e))^(1/2))*a^(1/2)*c^(1/2)/f-2*arcta
n(cos(f*x+e)*a^(1/2)*d^(1/2)/(a+a*sin(f*x+e))^(1/2)/(c+d*sin(f*x+e))^(1/2))*a^(1/2)*d^(1/2)/f

Rubi [A] (verified)

Time = 0.32 (sec) , antiderivative size = 123, normalized size of antiderivative = 1.00, number of steps used = 5, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.143, Rules used = {3028, 2854, 211, 3022, 212} \[ \int \csc (e+f x) \sqrt {a+a \sin (e+f x)} \sqrt {c+d \sin (e+f x)} \, dx=-\frac {2 \sqrt {a} \sqrt {d} \arctan \left (\frac {\sqrt {a} \sqrt {d} \cos (e+f x)}{\sqrt {a \sin (e+f x)+a} \sqrt {c+d \sin (e+f x)}}\right )}{f}-\frac {2 \sqrt {a} \sqrt {c} \text {arctanh}\left (\frac {\sqrt {a} \sqrt {c} \cos (e+f x)}{\sqrt {a \sin (e+f x)+a} \sqrt {c+d \sin (e+f x)}}\right )}{f} \]

[In]

Int[Csc[e + f*x]*Sqrt[a + a*Sin[e + f*x]]*Sqrt[c + d*Sin[e + f*x]],x]

[Out]

(-2*Sqrt[a]*Sqrt[d]*ArcTan[(Sqrt[a]*Sqrt[d]*Cos[e + f*x])/(Sqrt[a + a*Sin[e + f*x]]*Sqrt[c + d*Sin[e + f*x]])]
)/f - (2*Sqrt[a]*Sqrt[c]*ArcTanh[(Sqrt[a]*Sqrt[c]*Cos[e + f*x])/(Sqrt[a + a*Sin[e + f*x]]*Sqrt[c + d*Sin[e + f
*x]])])/f

Rule 211

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[a/b, 2]/a)*ArcTan[x/Rt[a/b, 2]], x] /; FreeQ[{a, b}, x]
&& PosQ[a/b]

Rule 212

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1/(Rt[a, 2]*Rt[-b, 2]))*ArcTanh[Rt[-b, 2]*(x/Rt[a, 2])], x]
/; FreeQ[{a, b}, x] && NegQ[a/b] && (GtQ[a, 0] || LtQ[b, 0])

Rule 2854

Int[Sqrt[(a_) + (b_.)*sin[(e_.) + (f_.)*(x_)]]/Sqrt[(c_.) + (d_.)*sin[(e_.) + (f_.)*(x_)]], x_Symbol] :> Dist[
-2*(b/f), Subst[Int[1/(b + d*x^2), x], x, b*(Cos[e + f*x]/(Sqrt[a + b*Sin[e + f*x]]*Sqrt[c + d*Sin[e + f*x]]))
], x] /; FreeQ[{a, b, c, d, e, f}, x] && NeQ[b*c - a*d, 0] && EqQ[a^2 - b^2, 0] && NeQ[c^2 - d^2, 0]

Rule 3022

Int[Sqrt[(a_) + (b_.)*sin[(e_.) + (f_.)*(x_)]]/(sin[(e_.) + (f_.)*(x_)]*Sqrt[(c_) + (d_.)*sin[(e_.) + (f_.)*(x
_)]]), x_Symbol] :> Dist[-2*(a/f), Subst[Int[1/(1 - a*c*x^2), x], x, Cos[e + f*x]/(Sqrt[a + b*Sin[e + f*x]]*Sq
rt[c + d*Sin[e + f*x]])], x] /; FreeQ[{a, b, c, d, e, f}, x] && NeQ[b*c - a*d, 0] && EqQ[a^2 - b^2, 0] && NeQ[
b*c + a*d, 0]

Rule 3028

Int[(Sqrt[(a_) + (b_.)*sin[(e_.) + (f_.)*(x_)]]*Sqrt[(c_) + (d_.)*sin[(e_.) + (f_.)*(x_)]])/sin[(e_.) + (f_.)*
(x_)], x_Symbol] :> Dist[d, Int[Sqrt[a + b*Sin[e + f*x]]/Sqrt[c + d*Sin[e + f*x]], x], x] + Dist[c, Int[Sqrt[a
+ b*Sin[e + f*x]]/(Sin[e + f*x]*Sqrt[c + d*Sin[e + f*x]]), x], x] /; FreeQ[{a, b, c, d, e, f}, x] && NeQ[b*c
- a*d, 0] && (NeQ[a^2 - b^2, 0] || NeQ[c^2 - d^2, 0])

Rubi steps \begin{align*} \text {integral}& = c \int \frac {\csc (e+f x) \sqrt {a+a \sin (e+f x)}}{\sqrt {c+d \sin (e+f x)}} \, dx+d \int \frac {\sqrt {a+a \sin (e+f x)}}{\sqrt {c+d \sin (e+f x)}} \, dx \\ & = -\frac {(2 a c) \text {Subst}\left (\int \frac {1}{1-a c x^2} \, dx,x,\frac {\cos (e+f x)}{\sqrt {a+a \sin (e+f x)} \sqrt {c+d \sin (e+f x)}}\right )}{f}-\frac {(2 a d) \text {Subst}\left (\int \frac {1}{a+d x^2} \, dx,x,\frac {a \cos (e+f x)}{\sqrt {a+a \sin (e+f x)} \sqrt {c+d \sin (e+f x)}}\right )}{f} \\ & = -\frac {2 \sqrt {a} \sqrt {d} \arctan \left (\frac {\sqrt {a} \sqrt {d} \cos (e+f x)}{\sqrt {a+a \sin (e+f x)} \sqrt {c+d \sin (e+f x)}}\right )}{f}-\frac {2 \sqrt {a} \sqrt {c} \text {arctanh}\left (\frac {\sqrt {a} \sqrt {c} \cos (e+f x)}{\sqrt {a+a \sin (e+f x)} \sqrt {c+d \sin (e+f x)}}\right )}{f} \\ \end{align*}

Mathematica [C] (warning: unable to verify)

Result contains complex when optimal does not.

Time = 2.90 (sec) , antiderivative size = 567, normalized size of antiderivative = 4.61 \[ \int \csc (e+f x) \sqrt {a+a \sin (e+f x)} \sqrt {c+d \sin (e+f x)} \, dx=-\frac {\left (\sqrt {c} \log \left (\frac {\left (\frac {1}{2}+\frac {i}{2}\right ) e^{-\frac {i e}{2}} \left (-\sqrt {2} c \left (-1+e^{i (e+f x)}\right )-i \sqrt {2} d \left (1+e^{i (e+f x)}\right )+2 i \sqrt {c} \sqrt {2 c e^{i (e+f x)}-i d \left (-1+e^{2 i (e+f x)}\right )}\right ) f}{c^{3/2} \left (1+e^{i (e+f x)}\right )}\right )+\sqrt {c} \log \left (\frac {\left (\frac {1}{2}+\frac {i}{2}\right ) e^{-\frac {i e}{2}} \left (-i \sqrt {2} d \left (-1+e^{i (e+f x)}\right )+\sqrt {2} c \left (1+e^{i (e+f x)}\right )+2 \sqrt {c} \sqrt {2 c e^{i (e+f x)}-i d \left (-1+e^{2 i (e+f x)}\right )}\right ) f}{c^{3/2} \left (-1+e^{i (e+f x)}\right )}\right )-i \sqrt {d} \left (\log \left (\frac {2 e^{-\frac {1}{2} i (e+2 f x)} \left ((-1)^{3/4} d+\sqrt [4]{-1} c e^{i (e+f x)}+i \sqrt {d} \sqrt {2 c e^{i (e+f x)}-i d \left (-1+e^{2 i (e+f x)}\right )}\right ) f}{d^{3/2}}\right )-\log \left (\frac {(1+i) \sqrt {2} \left (c-i d \cos (e+f x)+d \sin (e+f x)+(1-i) \sqrt {d} \sqrt {(\cos (e+f x)+i \sin (e+f x)) (c+d \sin (e+f x))}\right )}{\sqrt {d}}\right )\right )\right ) \left (\cos \left (\frac {1}{2} (e+f x)\right )+i \sin \left (\frac {1}{2} (e+f x)\right )\right ) \sqrt {a (1+\sin (e+f x))} \sqrt {c+d \sin (e+f x)}}{f \left (\cos \left (\frac {1}{2} (e+f x)\right )+\sin \left (\frac {1}{2} (e+f x)\right )\right ) \sqrt {(\cos (e+f x)+i \sin (e+f x)) (c+d \sin (e+f x))}} \]

[In]

Integrate[Csc[e + f*x]*Sqrt[a + a*Sin[e + f*x]]*Sqrt[c + d*Sin[e + f*x]],x]

[Out]

-(((Sqrt[c]*Log[((1/2 + I/2)*(-(Sqrt[2]*c*(-1 + E^(I*(e + f*x)))) - I*Sqrt[2]*d*(1 + E^(I*(e + f*x))) + (2*I)*
Sqrt[c]*Sqrt[2*c*E^(I*(e + f*x)) - I*d*(-1 + E^((2*I)*(e + f*x)))])*f)/(c^(3/2)*E^((I/2)*e)*(1 + E^(I*(e + f*x
))))] + Sqrt[c]*Log[((1/2 + I/2)*((-I)*Sqrt[2]*d*(-1 + E^(I*(e + f*x))) + Sqrt[2]*c*(1 + E^(I*(e + f*x))) + 2*
Sqrt[c]*Sqrt[2*c*E^(I*(e + f*x)) - I*d*(-1 + E^((2*I)*(e + f*x)))])*f)/(c^(3/2)*E^((I/2)*e)*(-1 + E^(I*(e + f*
x))))] - I*Sqrt[d]*(Log[(2*((-1)^(3/4)*d + (-1)^(1/4)*c*E^(I*(e + f*x)) + I*Sqrt[d]*Sqrt[2*c*E^(I*(e + f*x)) -
I*d*(-1 + E^((2*I)*(e + f*x)))])*f)/(d^(3/2)*E^((I/2)*(e + 2*f*x)))] - Log[((1 + I)*Sqrt[2]*(c - I*d*Cos[e +
f*x] + d*Sin[e + f*x] + (1 - I)*Sqrt[d]*Sqrt[(Cos[e + f*x] + I*Sin[e + f*x])*(c + d*Sin[e + f*x])]))/Sqrt[d]])
)*(Cos[(e + f*x)/2] + I*Sin[(e + f*x)/2])*Sqrt[a*(1 + Sin[e + f*x])]*Sqrt[c + d*Sin[e + f*x]])/(f*(Cos[(e + f*
x)/2] + Sin[(e + f*x)/2])*Sqrt[(Cos[e + f*x] + I*Sin[e + f*x])*(c + d*Sin[e + f*x])]))

Maple [B] (warning: unable to verify)

Leaf count of result is larger than twice the leaf count of optimal. \(5065\) vs. \(2(99)=198\).

Time = 0.84 (sec) , antiderivative size = 5066, normalized size of antiderivative = 41.19

\[\text {output too large to display}\]

[In]

int((a+a*sin(f*x+e))^(1/2)*(c+d*sin(f*x+e))^(1/2)/sin(f*x+e),x)

[Out]

result too large to display

Fricas [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 333 vs. \(2 (99) = 198\).

Time = 1.02 (sec) , antiderivative size = 3539, normalized size of antiderivative = 28.77 \[ \int \csc (e+f x) \sqrt {a+a \sin (e+f x)} \sqrt {c+d \sin (e+f x)} \, dx=\text {Too large to display} \]

[In]

integrate((a+a*sin(f*x+e))^(1/2)*(c+d*sin(f*x+e))^(1/2)/sin(f*x+e),x, algorithm="fricas")

[Out]

[1/4*(sqrt(a*c)*log(((a*c^4 - 28*a*c^3*d + 70*a*c^2*d^2 - 28*a*c*d^3 + a*d^4)*cos(f*x + e)^5 + a*c^4 + 4*a*c^3
*d + 6*a*c^2*d^2 + 4*a*c*d^3 + a*d^4 - (31*a*c^4 - 196*a*c^3*d + 154*a*c^2*d^2 - 4*a*c*d^3 - a*d^4)*cos(f*x +
e)^4 - 2*(81*a*c^4 - 252*a*c^3*d + 150*a*c^2*d^2 - 28*a*c*d^3 + a*d^4)*cos(f*x + e)^3 + 2*(79*a*c^4 - 100*a*c^
3*d + 74*a*c^2*d^2 - 4*a*c*d^3 - a*d^4)*cos(f*x + e)^2 - 8*((c^3 - 7*c^2*d + 7*c*d^2 - d^3)*cos(f*x + e)^4 - 2
*(5*c^3 - 14*c^2*d + 5*c*d^2)*cos(f*x + e)^3 + 51*c^3 - 59*c^2*d + 17*c*d^2 - d^3 - 2*(18*c^3 - 33*c^2*d + 12*
c*d^2 - d^3)*cos(f*x + e)^2 + 2*(13*c^3 - 14*c^2*d + 5*c*d^2)*cos(f*x + e) + ((c^3 - 7*c^2*d + 7*c*d^2 - d^3)*
cos(f*x + e)^3 - 51*c^3 + 59*c^2*d - 17*c*d^2 + d^3 + (11*c^3 - 35*c^2*d + 17*c*d^2 - d^3)*cos(f*x + e)^2 - (2
5*c^3 - 31*c^2*d + 7*c*d^2 - d^3)*cos(f*x + e))*sin(f*x + e))*sqrt(a*c)*sqrt(a*sin(f*x + e) + a)*sqrt(d*sin(f*
x + e) + c) + (289*a*c^4 - 476*a*c^3*d + 230*a*c^2*d^2 - 28*a*c*d^3 + a*d^4)*cos(f*x + e) + (a*c^4 + 4*a*c^3*d
+ 6*a*c^2*d^2 + 4*a*c*d^3 + a*d^4 + (a*c^4 - 28*a*c^3*d + 70*a*c^2*d^2 - 28*a*c*d^3 + a*d^4)*cos(f*x + e)^4 +
32*(a*c^4 - 7*a*c^3*d + 7*a*c^2*d^2 - a*c*d^3)*cos(f*x + e)^3 - 2*(65*a*c^4 - 140*a*c^3*d + 38*a*c^2*d^2 - 12
*a*c*d^3 + a*d^4)*cos(f*x + e)^2 - 32*(9*a*c^4 - 15*a*c^3*d + 7*a*c^2*d^2 - a*c*d^3)*cos(f*x + e))*sin(f*x + e
))/(cos(f*x + e)^5 + cos(f*x + e)^4 - 2*cos(f*x + e)^3 - 2*cos(f*x + e)^2 + (cos(f*x + e)^4 - 2*cos(f*x + e)^2
+ 1)*sin(f*x + e) + cos(f*x + e) + 1)) + sqrt(-a*d)*log((128*a*d^4*cos(f*x + e)^5 + a*c^4 + 4*a*c^3*d + 6*a*c
^2*d^2 + 4*a*c*d^3 + a*d^4 + 128*(2*a*c*d^3 - a*d^4)*cos(f*x + e)^4 - 32*(5*a*c^2*d^2 - 14*a*c*d^3 + 13*a*d^4)
*cos(f*x + e)^3 - 32*(a*c^3*d - 2*a*c^2*d^2 + 9*a*c*d^3 - 4*a*d^4)*cos(f*x + e)^2 - 8*(16*d^3*cos(f*x + e)^4 +
24*(c*d^2 - d^3)*cos(f*x + e)^3 - c^3 + 17*c^2*d - 59*c*d^2 + 51*d^3 - 2*(5*c^2*d - 26*c*d^2 + 33*d^3)*cos(f*
x + e)^2 - (c^3 - 7*c^2*d + 31*c*d^2 - 25*d^3)*cos(f*x + e) + (16*d^3*cos(f*x + e)^3 + c^3 - 17*c^2*d + 59*c*d
^2 - 51*d^3 - 8*(3*c*d^2 - 5*d^3)*cos(f*x + e)^2 - 2*(5*c^2*d - 14*c*d^2 + 13*d^3)*cos(f*x + e))*sin(f*x + e))
*sqrt(-a*d)*sqrt(a*sin(f*x + e) + a)*sqrt(d*sin(f*x + e) + c) + (a*c^4 - 28*a*c^3*d + 230*a*c^2*d^2 - 476*a*c*
d^3 + 289*a*d^4)*cos(f*x + e) + (128*a*d^4*cos(f*x + e)^4 + a*c^4 + 4*a*c^3*d + 6*a*c^2*d^2 + 4*a*c*d^3 + a*d^
4 - 256*(a*c*d^3 - a*d^4)*cos(f*x + e)^3 - 32*(5*a*c^2*d^2 - 6*a*c*d^3 + 5*a*d^4)*cos(f*x + e)^2 + 32*(a*c^3*d
- 7*a*c^2*d^2 + 15*a*c*d^3 - 9*a*d^4)*cos(f*x + e))*sin(f*x + e))/(cos(f*x + e) + sin(f*x + e) + 1)))/f, 1/4*
(2*sqrt(-a*c)*arctan(-1/4*((c^2 - 6*c*d + d^2)*cos(f*x + e)^2 - 9*c^2 + 6*c*d - d^2 + 8*(c^2 - c*d)*sin(f*x +
e))*sqrt(-a*c)*sqrt(a*sin(f*x + e) + a)*sqrt(d*sin(f*x + e) + c)/((a*c^2*d - a*c*d^2)*cos(f*x + e)^3 - (a*c^3
- 3*a*c^2*d)*cos(f*x + e)*sin(f*x + e) + (2*a*c^3 - a*c^2*d + a*c*d^2)*cos(f*x + e))) + sqrt(-a*d)*log((128*a*
d^4*cos(f*x + e)^5 + a*c^4 + 4*a*c^3*d + 6*a*c^2*d^2 + 4*a*c*d^3 + a*d^4 + 128*(2*a*c*d^3 - a*d^4)*cos(f*x + e
)^4 - 32*(5*a*c^2*d^2 - 14*a*c*d^3 + 13*a*d^4)*cos(f*x + e)^3 - 32*(a*c^3*d - 2*a*c^2*d^2 + 9*a*c*d^3 - 4*a*d^
4)*cos(f*x + e)^2 - 8*(16*d^3*cos(f*x + e)^4 + 24*(c*d^2 - d^3)*cos(f*x + e)^3 - c^3 + 17*c^2*d - 59*c*d^2 + 5
1*d^3 - 2*(5*c^2*d - 26*c*d^2 + 33*d^3)*cos(f*x + e)^2 - (c^3 - 7*c^2*d + 31*c*d^2 - 25*d^3)*cos(f*x + e) + (1
6*d^3*cos(f*x + e)^3 + c^3 - 17*c^2*d + 59*c*d^2 - 51*d^3 - 8*(3*c*d^2 - 5*d^3)*cos(f*x + e)^2 - 2*(5*c^2*d -
14*c*d^2 + 13*d^3)*cos(f*x + e))*sin(f*x + e))*sqrt(-a*d)*sqrt(a*sin(f*x + e) + a)*sqrt(d*sin(f*x + e) + c) +
(a*c^4 - 28*a*c^3*d + 230*a*c^2*d^2 - 476*a*c*d^3 + 289*a*d^4)*cos(f*x + e) + (128*a*d^4*cos(f*x + e)^4 + a*c^
4 + 4*a*c^3*d + 6*a*c^2*d^2 + 4*a*c*d^3 + a*d^4 - 256*(a*c*d^3 - a*d^4)*cos(f*x + e)^3 - 32*(5*a*c^2*d^2 - 6*a
*c*d^3 + 5*a*d^4)*cos(f*x + e)^2 + 32*(a*c^3*d - 7*a*c^2*d^2 + 15*a*c*d^3 - 9*a*d^4)*cos(f*x + e))*sin(f*x + e
))/(cos(f*x + e) + sin(f*x + e) + 1)))/f, 1/4*(2*sqrt(a*d)*arctan(1/4*(8*d^2*cos(f*x + e)^2 - c^2 + 6*c*d - 9*
d^2 - 8*(c*d - d^2)*sin(f*x + e))*sqrt(a*d)*sqrt(a*sin(f*x + e) + a)*sqrt(d*sin(f*x + e) + c)/(2*a*d^3*cos(f*x
+ e)^3 - (3*a*c*d^2 - a*d^3)*cos(f*x + e)*sin(f*x + e) - (a*c^2*d - a*c*d^2 + 2*a*d^3)*cos(f*x + e))) + sqrt(
a*c)*log(((a*c^4 - 28*a*c^3*d + 70*a*c^2*d^2 - 28*a*c*d^3 + a*d^4)*cos(f*x + e)^5 + a*c^4 + 4*a*c^3*d + 6*a*c^
2*d^2 + 4*a*c*d^3 + a*d^4 - (31*a*c^4 - 196*a*c^3*d + 154*a*c^2*d^2 - 4*a*c*d^3 - a*d^4)*cos(f*x + e)^4 - 2*(8
1*a*c^4 - 252*a*c^3*d + 150*a*c^2*d^2 - 28*a*c*d^3 + a*d^4)*cos(f*x + e)^3 + 2*(79*a*c^4 - 100*a*c^3*d + 74*a*
c^2*d^2 - 4*a*c*d^3 - a*d^4)*cos(f*x + e)^2 - 8*((c^3 - 7*c^2*d + 7*c*d^2 - d^3)*cos(f*x + e)^4 - 2*(5*c^3 - 1
4*c^2*d + 5*c*d^2)*cos(f*x + e)^3 + 51*c^3 - 59*c^2*d + 17*c*d^2 - d^3 - 2*(18*c^3 - 33*c^2*d + 12*c*d^2 - d^3
)*cos(f*x + e)^2 + 2*(13*c^3 - 14*c^2*d + 5*c*d^2)*cos(f*x + e) + ((c^3 - 7*c^2*d + 7*c*d^2 - d^3)*cos(f*x + e
)^3 - 51*c^3 + 59*c^2*d - 17*c*d^2 + d^3 + (11*c^3 - 35*c^2*d + 17*c*d^2 - d^3)*cos(f*x + e)^2 - (25*c^3 - 31*
c^2*d + 7*c*d^2 - d^3)*cos(f*x + e))*sin(f*x + e))*sqrt(a*c)*sqrt(a*sin(f*x + e) + a)*sqrt(d*sin(f*x + e) + c)
+ (289*a*c^4 - 476*a*c^3*d + 230*a*c^2*d^2 - 28*a*c*d^3 + a*d^4)*cos(f*x + e) + (a*c^4 + 4*a*c^3*d + 6*a*c^2*
d^2 + 4*a*c*d^3 + a*d^4 + (a*c^4 - 28*a*c^3*d + 70*a*c^2*d^2 - 28*a*c*d^3 + a*d^4)*cos(f*x + e)^4 + 32*(a*c^4
- 7*a*c^3*d + 7*a*c^2*d^2 - a*c*d^3)*cos(f*x + e)^3 - 2*(65*a*c^4 - 140*a*c^3*d + 38*a*c^2*d^2 - 12*a*c*d^3 +
a*d^4)*cos(f*x + e)^2 - 32*(9*a*c^4 - 15*a*c^3*d + 7*a*c^2*d^2 - a*c*d^3)*cos(f*x + e))*sin(f*x + e))/(cos(f*x
+ e)^5 + cos(f*x + e)^4 - 2*cos(f*x + e)^3 - 2*cos(f*x + e)^2 + (cos(f*x + e)^4 - 2*cos(f*x + e)^2 + 1)*sin(f
*x + e) + cos(f*x + e) + 1)))/f, 1/2*(sqrt(-a*c)*arctan(-1/4*((c^2 - 6*c*d + d^2)*cos(f*x + e)^2 - 9*c^2 + 6*c
*d - d^2 + 8*(c^2 - c*d)*sin(f*x + e))*sqrt(-a*c)*sqrt(a*sin(f*x + e) + a)*sqrt(d*sin(f*x + e) + c)/((a*c^2*d
- a*c*d^2)*cos(f*x + e)^3 - (a*c^3 - 3*a*c^2*d)*cos(f*x + e)*sin(f*x + e) + (2*a*c^3 - a*c^2*d + a*c*d^2)*cos(
f*x + e))) + sqrt(a*d)*arctan(1/4*(8*d^2*cos(f*x + e)^2 - c^2 + 6*c*d - 9*d^2 - 8*(c*d - d^2)*sin(f*x + e))*sq
rt(a*d)*sqrt(a*sin(f*x + e) + a)*sqrt(d*sin(f*x + e) + c)/(2*a*d^3*cos(f*x + e)^3 - (3*a*c*d^2 - a*d^3)*cos(f*
x + e)*sin(f*x + e) - (a*c^2*d - a*c*d^2 + 2*a*d^3)*cos(f*x + e))))/f]

Sympy [F]

\[ \int \csc (e+f x) \sqrt {a+a \sin (e+f x)} \sqrt {c+d \sin (e+f x)} \, dx=\int \frac {\sqrt {a \left (\sin {\left (e + f x \right )} + 1\right )} \sqrt {c + d \sin {\left (e + f x \right )}}}{\sin {\left (e + f x \right )}}\, dx \]

[In]

integrate((a+a*sin(f*x+e))**(1/2)*(c+d*sin(f*x+e))**(1/2)/sin(f*x+e),x)

[Out]

Integral(sqrt(a*(sin(e + f*x) + 1))*sqrt(c + d*sin(e + f*x))/sin(e + f*x), x)

Maxima [F]

\[ \int \csc (e+f x) \sqrt {a+a \sin (e+f x)} \sqrt {c+d \sin (e+f x)} \, dx=\int { \frac {\sqrt {a \sin \left (f x + e\right ) + a} \sqrt {d \sin \left (f x + e\right ) + c}}{\sin \left (f x + e\right )} \,d x } \]

[In]

integrate((a+a*sin(f*x+e))^(1/2)*(c+d*sin(f*x+e))^(1/2)/sin(f*x+e),x, algorithm="maxima")

[Out]

integrate(sqrt(a*sin(f*x + e) + a)*sqrt(d*sin(f*x + e) + c)/sin(f*x + e), x)

Giac [F(-1)]

Timed out. \[ \int \csc (e+f x) \sqrt {a+a \sin (e+f x)} \sqrt {c+d \sin (e+f x)} \, dx=\text {Timed out} \]

[In]

integrate((a+a*sin(f*x+e))^(1/2)*(c+d*sin(f*x+e))^(1/2)/sin(f*x+e),x, algorithm="giac")

[Out]

Timed out

Mupad [F(-1)]

Timed out. \[ \int \csc (e+f x) \sqrt {a+a \sin (e+f x)} \sqrt {c+d \sin (e+f x)} \, dx=\int \frac {\sqrt {a+a\,\sin \left (e+f\,x\right )}\,\sqrt {c+d\,\sin \left (e+f\,x\right )}}{\sin \left (e+f\,x\right )} \,d x \]

[In]

int(((a + a*sin(e + f*x))^(1/2)*(c + d*sin(e + f*x))^(1/2))/sin(e + f*x),x)

[Out]

int(((a + a*sin(e + f*x))^(1/2)*(c + d*sin(e + f*x))^(1/2))/sin(e + f*x), x)